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We propose a new method to compute the running coupling constant of gauge theories on the 
lattice. We first give the definition of the running coupling in the new scheme using the Wilson 
loops in a finite volume, and explain how the running of the coupling constant is extracted from 
the measurement of the volume dependence. The perturbative calculation of the renormalization 
constant to define the scheme is also given at the leading order. As a benchmark test of the new 
scheme we apply the method the case of the quenched QCD. We show the preliminary result from 
our numerical simulations which are carried out with plaquette gauge action for various lattice 
sizes and bare lattice couplings. With techniques to improve the statistical accuracy, we show that 
we can determine the non-perturbative running of the coupling constant in a wide range of the 
energy scale with relatively small number of gauge configurations in our scheme. We compare our 
lattice data of the running coupling constant with perturbative renormalization group evolution at 
one- and two-loop order, and confirm the consistency between them at high energy. 
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1. Introduction 

The properties of a vectorial gauge theory which has (near) conformal infrared fixed point 
(IRFP) are of fundamental importance. In addition to its intrinsic field-theoretic interest, this ap- 
proximate conformal, or "walking", behavior is an essential ingredient of modern technicolor mod- 
els of dynamical electroweak symmetry breaking [p, providing requisite enhancement of Standard 
Model fermion masses [g, Q, g, 0]. This walking behavior could also give a remedy for the 
problem of large correction (often denoted as S [|| ^]) to the Z boson propagator |[T^, 11,12]. Our 



ultimate objective is to investigate the nature of physics in such a near conformal gauge theory. 
However, we first have to address the following question: is there any theory which actually has a 
conformal nature in a certain energy region? 

An SU (N) gauge theory with a large number of massless fermions has been known as a 



promising candidate for a theory which has the nontrivial IRFP [13]. For example, in the case 



of SU (3) gauge theory, when the number of massless fermions, Nf, takes values in the range 

4^(33/2-^/) 
19(iV/-153/19) ' 



8 < Nf < 16, the two-loop running coupling approaches to a finite value, a* = < i n 



the infrared energy region due to the existence of the IRFP. When Nf is close to 16 (which is the 
largest Nf we can take when we restrict our study to asymptotically free theories), even though the 
value of a* is corrected by higher order terms in the perturbative expansion, the existence of the 
IRFP is valid beyond two-loop approximation since the coupling is small enough. However, when 
we consider a range of Nf where the above two-loop estimation of a* becomes large, it is quite 
non-trivial whether this IRFP really exists beyond perturbation theory. 

In principle, lattice simulations should provide a way to determine whether an SU (N) gauge 
theory with a certain number of massless fermions has conformal fixed point or not. Several pio- 



neering works Q14Q found the chiral phase transition at certain critical flavor, which can be consid- 
ered as an indication of the appearance of the IRFP. However, groups that have studied this have 
not reached consensus on the critical number of flavor for the chiral phase transition. Recently, Ap- 



pelquist et al. Ql 5[ | directly calculated the running coupling of the SU (3) gauge theory with 8 and 
12 flavors by using the Schrodinger functional (SF) method on the lattice, and showed an evidence 
for the existence of the IRFP for Nf = 12 while the running coupling for Nf = 8 is more like QCD 
with three flavors. 

Here, we propose a new scheme for the calculation of the running coupling on the lattice. 
We extract the renormalized coupling from the Wilson loops in a finite volume, and determine 
the running of the coupling constant from the measurement of the volume dependence by using 
the step scaling procedure. This new scheme can serve as an alternative computational method in 
general gauge theories. 

In the context of exploring the conformal window, using a different scheme to calculate the 
coupling constant is well-motivated by itself to verify the real existence of the fixed point. Technical 
advantage of our scheme is that it is free from 0(a) error in contrast to the SF scheme, in which it 
is practically difficult to remove all the 0(a) discretization effects. 

In Section 2, we define the new scheme for the determination of the running coupling. Then 
we review the step-scaling procedure to show how to actually obtain the evolution of the running 
coupling from lattice simulations in Section 3. Sections 4 and 5 are devoted to the explanation 
of the details of our numerical simulations, and in Section 6, we show the result of the quenched 



2 



A new method of calculating the running coupling constant 



study as a test to see the effectiveness of our new scheme for the calculation of the running cou- 
pling. Discussion on the numerical results are also given in Section 6. Section 7 summarizes our 
conclusions. 



2. Wilson Loop Scheme 



In this section, we give the definition of the new renormalization scheme, the "Wilson loop 
scheme", for the running coupling, and show how to calculate it on the lattice. Let us start with 
general features in the renormalization of the coupling constant. Consider an amplitude A whose 
tree-level contribution is proportional to gf } (where go is the bare coupling constant): 



(2.1) 



Here, k is a certain coefficient which is a function of all the parameters of the theory except go- 
Then, we denote the ratio of the fully non-perturbative value of the amplitude A to its tree-level 
value as Z(/i): 

(2.2) 



A NP (/i) = Z( J u)A tl 



where jj. is the scale at which the amplitude A is defined. By using Eq. < \2.l\ ), the right hand side of 
the above equation can be rewritten as Z(/l)gQ k, and the combination Z(}i,)g^ can be identified as 
the renormalized coupling at the scale pL. So the renormalized coupling, g(ju), can be expressed as 
follows: 



(2.3) 



Thus, any amplitude with a tree-level value proportional to g^ can be used to define the renormal- 
ized coupling. Here, we use the following quantity: 



dRdT 



ln(W(R,T;L ,T )) 



(2.4) 



T=R 



where W(R, T) is the Wilson loop. The definition of the Wilson loop is graphically shown in Fig. |Tj. 
In this figure, To, Lq, and T , R represent the size of the box and the Wilson loop in the temporal 
and spatial directions, respectively, and a is the lattice spacing. From now on, for simplicity, we 
consider the case of 7b = Lq. At tree level in the perturbative expansion, this quantity actually is 
proportional to gjj, i.e., 

92 : ln(W(/?,r;L )) tree 



-R z 



dRdT 



(2.5) 



T=R 



where k, in the case of periodic boundary condition for example, is 



-R 1 



dRdT (2k) 4 
+ zero mode contribution. 



I 
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T=R 



(2.6) 



Here, (no, «i ,12,13) represents integer four-vector to define the momentum. "Zero mode contribu- 
tion" in the above equation is coming from the existence of so-called "toron" contributions which 
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Figure 1: Wilson loop defined on the latticized space-time box. 7b, Lq and T, R represent the size of the box 
and the Wilson loop in the temporal and spatial directions, respectively, a is the lattice spacing. 
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Figure 2: Values of k for several values of R/Lq and Lo/a (colored squares whose Lo/ais indicated by the 
numbers with the same color). The value of k in the continuum limit is also shown as a solid curve. 



originate from zero-mode configurations degenerate with the vacuum on the periodic torus. This 



contribution is calculated in Ref. []16[], and we use the result from that paper. After evaluating the 
summation in Eq. (2.6) , one can find that k only depends on the value of R/Lq. The value of k as 
a function of R/Lq in the continuum limit is shown in Fig. ||. We also did similar calculations of k 
in the case of discrete space-time, and plotted them for several values of Lo/a and R/Lq. Note that 
the continuum limit actually exists {i.e., k is finite in the limit of Lq/o — > °°) and that the conver- 
gence to continuum value is faster for larger values of R/Lq. Once the value of k is obtained, the 



'Detailed calculation of the factor k can be found in Jl7|]. 
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renormalized coupling is denned according to Eq. (2.3): 

(/?,r;L )) NP 

(2.7) 



R\ -R 2 ^fln(W(R,T;L )Y 
g z [ Lq,— 



The numerator on the right hand side of Eq. ( |2.7| ) can be estimated from the Creutz ratio on the 
lattice, 

+ 1/2, f + l/2;Lo/a) = -In ~ f ^ f , (2.8) 

\W(R+\,T;L /a) W(R,T + l;L /a) J 

where f = T/a and R = R/a. The value of £ is evaluated by a Monte Carlo (MC) simulation. 
Then the renormalized coupling constant in the Wilson loop scheme can be written as: 

^flo,*±^,± ) =(R+l/2) 2 . X (R+l/2;L Q /a)/k, (2.9) 

V M) M) / 

where we used the shorthand notation ;£()? + 1/2,1" + l/2;Lo/a)|j$ = f = x(R+ l/2;Lo/a). Here, 
since depends on three different scales, namely, Lo, R, and a, we indicated that can be 
viewed as a function of Lo, (R + a/2) /Lq(= r), and a/Lo. We choose a specific value of r (r = 0.3, 
for example) and keep it fixed to that value throughout the analysis. Varying r means changing 
renormalization scheme. As for a/Lo, we extrapolate it to zero when taking the continuum limit. 
After fixing these two dimensionless parameters r and a/L , becomes a function of only one 
scale, Lo. In our scheme, Lo is identified as the scale at which the renormalized coupling is defined. 

One thing we should emphasize here is that the Creutz ratio is free from 0(a) discretization 
error, mainly because 0(a) -improvement of the heavy quark propagator is automatically achieved 
after the redefinition of the mass and the wavefunction [|l^]. Thus, our scheme explained here does 
not have any 0(a) systematic error as long as we use actions which do not have 0(a) error. 



3. Step scaling 

Here, we review the step-scaling procedure [jl9, 20] in the Wilson loop scheme, which enables 
us to evaluate the evolution of the running coupling for a large range of energy scale on the lattice. 
First, we choose a specific value for g^, g^ = g^. Then, for a fixed value of r, we find sets of 
parameters, Q3,Lo/a), which reproduce for several different values of Lo/a: 

{[pi 1 \(Lo/a)?),[^\(L /a)^),---}. (3.1) 

What we are doing here is tuning the value of /3 in such a way that the physical volume Lo is fixed 
for different values of La/ a. Let us call this fixed physical volume for the starting point as Lo, 
i.e., gw(Lo) = g„ ■ Next thing to do is to vary the physical volume from Lo to sLo, which gives 
the evolution of the running coupling from the energy scale 1/Lo to I/sLq. Here s is the scaling 
factor. This step can be achieved by changing the lattice size from (L /a)W to s(L /a) (1) with 
each value of jS^ unchanged. Values of g\, calculated with these new parameter sets should be 
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considered as the coupling at the energy scale 1 /sLq up to discretization error, so the extrapolation 
to the continuum limit can be taken by using those data. 

1 a 



j(a) 



(3.2) 



jL Q sL Q/ 

The resultant value of coupling, g 2 R , should be considered as the renormalized coupling at the energy 
scale 1/sLq. This is the way to obtain a single discrete step of evolution of the running coupling 
with the scaling factor s. 

Next, we find new parameter sets of (f5^ 2 \(Lo/a)^) which reproduce the value of ^{sLq) 
obtained in the previous step. Here, we find these parameter sets in such a way that the new lattice 
size Lo/cffl is equal to the original small lattice {Lq/o)^ . From here, we can repeat exactly the 
same procedure as that of the first step scaling explained in the previous paragraph: we calculate 
g„ with the parameter sets (f$( 2 \s(Lo/a)W). By iterating this procedure, say, n times, we obtain 
the evolution of the running coupling from the energy scale 1/Lq to 1/ (s n Lo). 

In our exploratory study, we slightly modify the step scaling procedure. Instead of tuning 



j3 to keep g^(Lo) fixed, we use the results given by Alpha collaboration [21], namely, we use the 
parameter sets which are tuned to keep the renormalized coupling in the SF scheme, g| F (Lo), fixed. 



We also use the result given in ||22|] to obtain additional parameters for the study of lower energy 
region. 

4. Simulation parameters 

We consider a four-dimensional Euclidean lattice Lq. In this work, we use the standard Wilson 
plaquette gauge action. Gauge configurations are generated by the pseudo-heatbath algorithm and 
over-relaxation, mixed in the ratio of 1:5. We call the combination of one pseudo-heatbath update 
sweep followed by five over-relaxation sweeps "an iteration". In order to ensure decorrelation, we 
save gauge configurations separately by 1000 iterations and perform measurements on them. The 
number of gauge configurations we generated is 100 at each set of (j3,Lo/a). We use both periodic 
and twisted boundary conditions. However, in this article we only report the case of periodic 
boundary condition. 

We have five parameter sets (Set 1- Set 5) of (j8,Lo/a) corresponding to the fixed physical box 
size Lq as in Table [j]. 

The values of /3 for Set 1-4 are tuned by fixing Lq using the renormalized coupling in SF 



scheme, g SF , as inputs. (See, Table 6 in [21].) Also, these parameters are tuned in such a way 
that the value of g 2 SF obtained from the first column {i.e. s = 1) of Set 2 (or 3, 4) coincide with 
that obtained from the second column {i.e. s = 2) of Set 1 (or 2, 3) up to systematic and statistical 
errors. Thus, if we call the energy scale defined by the s = 1 of Set 1 as 1 /Lq, from the simulations 
with parameters in Set 1-4, we obtain the evolution of the running coupling from 1 /Lq to 1/ (2 4 Lq). 
We also added an extra set (Set 5) of parameters to obtain further step scalings toward low energy 
region. Parameters in Set 5 are tuned by using the Sommer scale, ro, as inputs. (See Eq. (2.18) in 



Note that even though we use existing parameter sets which are tuned by using g 2 SF or ro as an 
input here, to make our study self-contained, we should tune parameters by using g 2 , itself. Such 
study is underway, and will be presented in our future publication [|l7j] . 
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Set 1 






Set 2 






Set 3 






Set 4 






P 


Lo/a 


Lo/a 


P 


Lo/a 


L Q /a 


P 


L /a 


Lo/a 


P 


L /a 


Lo/a 




(5=1) 


(, = 2) 




(5=1) 


(5 = 2) 




( s =l) 


(s = 2) 




(5=1) 


(s = 2) 


8.2500 


(8) 


16 


7.6547 


(8) 


16 


7.0197 


(8) 


16 


6.4527 


(8) 


16 


8.4677 


(10) 


20 


7.8500 


(10) 


20 


7.2098 


(10) 


20 


6.6629 


(10) 


20 


8.5985 


12 


24 


7.9993 


12 


24 


7.3551 


12 


24 


6.7750 


12 


24 


8.7289 


14 




8.1352 


14 




7.4986 


14 




6.9169 


14 




8.8323 


16 




8.2415 


16 




7.0203 


16 




7.6101 


16 





Set 5 








P 


Lo/a 


Lo/a 


Lo/a 




(,= 1) 


(5=1.5) 


(5 = 2) 


6.1274 


(8) 


12 


16 


6.2647 


(10) 




20 


6.3831 


12 


18 


24 


6.4841 


14 






6.5700 


16 


24 





Table 1: The parameter sets of j3 and Lo/a used for measurements. Each of the first columns in Sets 1-4 
gives the constant SF coupling, and the first one in Set 5 gives the constant Sommer scale. The parameter 
sets of Lo/a — 8 and 10 (denoted by parentheses) are used only for the reference to set the scale, but not for 
measurements. 



5. Simulation details 



Now, we calculate the renormalized coupling constant in the Wilson loop scheme defined by 



Eq. (2.9) for a fixed value of r: 

gl( L o,a/L 



(R 2 + \/2) 2 - X (R+\/2;Lo/a)/k\ 



(5.1) 



First, to reduce the statistical error, we use the APE smearing Q23p of link variables defined by 
the following equation; 



U 



(n+l) 



Pr ojsu(3) 



An) 



1. 



r(»)rr(") 



r(»)t 



77 v > J- —Y^ ttv'Jttv'J 77 W 
u x,n i c ^ii^v u x,v u x+v,ii u x+ii, 



(5.2) 



where n and c denote a smearing level and a smearing parameter, respectively. The result does not 
depend on the value of c significantly, and we take c = 2.3 in this work. Note that we need to find 
a optimal set of r = R+ ^j 2 and the smearing level n by considering the following requirements. 
To control the discretization error, it is better to choose a larger r. For the purpose of reducing 
the statistical error, it is better to take a smaller r and higher n. Fig. |3| shows the smearing-level 
dependence of (R + 1 /2) 2 • % in the case of j3 = 6.3831 and Lo/a = 18 as an example. From this 
figure, we find the statistical error is notably reduced even at smearing level one. Furthermore, 
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in order to avoid over-smearing, n should be smaller than R/2. This condition gives the lower 
bound, Lo/a > (4n+ l)/2. We summarize the bound from this requirement in Tabled We actually 
see, for example in the case of Lo/a = 18 (see Fig. ||), that the data of (R+ 1/2) = 1.5, 2.5 in 
higher smearing level are not reliable because of over-smearing. By considering all of the above 
requirements, we find that (r,n) = (0.3, 1) is an optimal choice. 

(R+l/2) 2 - X 

3.0 
2.5 
2.0 
1 .5 
1 .0 



0.5 
0.0 

2 4 6 8 10 

£+1/2 

Figure 3: The values of (R + 1 /2) 2 • % with statistical error for several (R + 1 /2) in the case of j3 = 6.3831 
and Lo/a — 18. The black cross, the red square and the blue diamond denote the data with 0, 1 and 2 
smearing, respectively. 



n 


r = 0.25 


r = 0.30 


r = 0.35 


n=l 


L /a > 10 


Lo/a > 8.3 


Lo/a > 7.1 


n = 2 


L Q /a > 18 


Lo/a > 15 


Lo/a > 12.8 


71 = 3 


Lq /a > 26 


Lo/a > 21.6 


Lo/a > 18.5 



Table 2: The lower bound for Lo/a to avoid over smearing. 



From here, we fix the renormalization condition r = 0.3, then we have to calculate the value 
of (R + 1 /2) 2 ■ x for noninteger R. We interpolate the value of (R + 1 /2) 2 • % using a quadratic fit 
function: 

f(R+l/2)=co + c 1 {R + \/2)+c 2 {R + l/2) 2 . (5.3) 

We determine the fit ranges for each lattice size as in Table |3| We confirmed that the data can be 
well fitted by our fit function with these fit ranges for all parameter sets. An example is shown in 
Fig-B 

Finally, we take the continuum limit of the running coupling constant of each column in Ta- 
ble U We show the example of the extrapolation of the coupling constant in Fig. ||. The red line 
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L Q /a 


£+1/2 




Rmax 


12 


3.6 


2 


4 


14 


4.2 


2 


5 


16 


4.8 


3 


5 


18 


5.4 


4 


6 


20 


6.0 


4 


6 


24 


7.2 


5 


7 



Table 3: The fit range used to interpolate the value of (R + 1 /2) 2 • %. The column "R + 1 /2" is the value that 
corresponds to r = 0.3. 



(R+l/2) 2 -x 



0.8 



0.6 



0.2 



0.0 



£+1/2 



Figure 4: The plot for the interpolation of (R + 1 /2) 2 • £ in the case of j3 = 6.5700 and Lq /a = 16. The dia- 
monds show the resultant values of (R + 1 /2) 2 • % for integer R with statistical error. The red line represents 
the function f(R + 1 /2) fitted to the data in the range of 3 < R < 5. The blue box indicates the interpolated 
value corresponding to r = 0.3. 



corresponds to the extrapolation of the first column of Set 1, and the blue line is that of the second 
column of Set 1. The values at the continuum limit are the renormalized coupling constants at 
scales 1/Lo and l/(2Lo), respectively. In performing the continuum extrapolation, we have used 
the fit function linear in (a/Lo) 2 : 



Co+Cl ' t) ■ (5 " 4) 



As explained in section 2, our Wilson loop scheme does not have 0(a) systematic error. However, 
in the present quenched QCD test, since we took the SF coupling, which has 0(a) error, as a 
constant input parameter, it is possible that an 0(a) error is introduced also in our data. Fig. || 



shows that the data are nicely fitted by the linear function, Eq. (5.4). This would indicates that 
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0(a) systematic error and possible higher order discretization error are small as we expected from 



perturbative analysis Q24Q in this weak coupling regime. 



2.2 



2.0 



2 1 .8 



1 .4 



1 .2 




(a/L Q ) 2 

Figure 5: The continuum limit of g\, in Set 1. The red and blue lines describe the continuum limit of the 
first and second columns of Set 1, respectively. 



6. Results and discussion 



We plot the numerical results for the energy scale dependence of the renormalized coupling 
in the Wilson loop scheme in Fig. ^. The perturbative renormalization group evolution at one-loop 
and two-loop are also displayed in the figure for comparison. Our MC data are consistent with the 
perturbative results in high energy region, while there is difference in low energy region because 
of the higher loop or the nonperturbative effect. 

Recall that the second column (s = 2) of Set 1 (or 2, 3, 4) and the first column (s = 1) of 
Set 2 (or 3, 4, 5) in Table [I] give the same value of the renormalized coupling in the SF scheme 
within the statistical error. From Fig. ||, we see that our results in the Wilson loop scheme also 
show that the step scaling procedure works well; steps are nicely connected as in the case of the 
SF couplings. It suggests that the systematic errors, such as the interpolation in R+ 1/2, 0(a) 
error in SF coupling and the extrapolation to the continuum limit are under control. To make our 
study self-contained, we should extend this scheme by taking g 2 „ which has no 0(a) systematic 
error, as an input physical parameter. Such study is underway, and will be presented in our future 
publication [|7p. 

Finally, we emphasize that there is an optimal choice of smearing level and r with which both 
the statistical and discretization errors are under control. The number of gauge configurations we 
used for the present study was only 100 for each (jS,Lo/a). From this small number of gauge 
configurations, we obtained data with a few (~ 10) percent statistical error in high (low) energy 
region. Thus, we could expect that this scheme is promising for the study of other theories. 
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Figure 6: The running coupling in the Wilson loop scheme. The horizontal axis is the ratio of the energy 
scale Lq/Lq, and the vertical axis is the renormalized coupling in the Wilson loop scheme. The square 
symbols indicate our lattice data of the running coupling constants. The right-most symbol corresponds to 
the renormalized coupling given by the first column of Set 1 at scale 1/Lq, the adjacent two overlapping 
symbols give the couplings from the second column of Set 1 and the first column of Set 2 at scale 1 / (2Lq), 
and similarly for the remaining data for lower energy scales. The perturbative running couplings at one-loop 
(dotted line) and two-loop (solid line) are also shown for comparison. 

7. Summary 

We proposed a new scheme for the determination of the running coupling on the lattice. Our 
method is based on the measurement of the finite volume dependence of the Wilson loop. Unlike 
the SF scheme, our method does not have any 0(a) discretization error, therefore the systematic er- 
ror arising from the extrapolation to the continuum limit is expected to be quite small. We showed 
results of preliminary numerical study for the quenched QCD as a feasibility test of our scheme, and 
confirmed that the method actually reproduced the step scaling of the coupling which is consistent 
with the two-loop running coupling at high energy. We also showed that the coupling calculated 
by this newly proposed scheme deviates from that with two-loop approximation below a certain 
energy scale. This deviation arises from the effects that are not captured by the two-loop approxi- 
mation. We have confirmed that our scheme works well for the calculation of the running coupling 
with relatively small number of gauge configurations, via demonstrating that the statistical error is 
under control by properly choosing the smearing level and r. We expect that this new method is 
quite useful to calculate the running coupling of the SU (N) gauge theory with a large number of 
dynamical fermions, which will be studied in our future work. 

The properties of a vectorial gauge theory as a function of the number of fermions is quite 
interesting. Through the analyses based on the Schwinger-Dyson and Bethe-Salpeter equations, it 
is shown that there is chiral symmetry restoration at certain critical number of fermions [||, [7p, and 
that the behavior of the meson spectrum, the S parameter, etc. , are quite different from those in 
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QCD with three flavors. [[II], [12|, ^5|, pffi ). In these analyses, two-loop running coupling was used 
as an approximation to the fully non-perturbative vertices. Therefore, the existence of the IRFP is 
an assumption and an input in these analyses. The confirmation of the existence of the IRFP by 
lattice simulations justifies this assumption, and motivates a fully non-perturbative study of near 
conformal gauge theories on the lattice. 
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